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In  this  paper  we  propose  a  new  way  for  Mossbauer  data  treatment  when  numerous  experimental  spectra 
are  recorded  in  operando  conditions  depending  on  a  perturbation.  In  our  example,  the  perturbation  is 
the  Li  amount  extracted  from  a  positive  electrode  material  LiFeo.75Mn0.25P04.  In  other  cases  perturbation 
could  be  the  recording  temperature,  the  pressure  or  kinetic  parameter  or  even  time  for  isothermal  exper¬ 
iments.  From  analysis  of  both  synchronous  and  asynchronous  2D-correlation  spectra,  we  can  focus  our 
attention  on  the  intensity  variations  at  some  specific  positions  deduced  from  cross-peaks.  This  analytical 
method  is  very  powerful  when  overlap  between  absorption  lines  is  observed.  This  is  typically  the  case 
when  Fe2+/Fe3+  contributions  are  simultaneously  present  in  a  Mossbauer  spectrum  at  lower  velocities. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  electrochemical  reactions  of  battery  materials  have  been 
extensively  studied  by  various  electrochemical  and  spectroscopic 
techniques  [1,2].  We  recently  determined  Lamb-Mossbauer  fac¬ 
tors  of  both  FeP04  and  LiFeP04  compounds  [3]  in  order  to  get 
precise  contributions  of  each  formed/reacted  species  during  elec¬ 
trochemical  process.  Mossbauer  is  a  powerful  technique  to  follow 
electrochemical  reactions  studied  by  ex  situ  [4,5],  in  situ  [6]  or 
operando  conditions  and  also  for  characterization  of  pristine  mate¬ 
rials.  In  some  case  it  is  possible  to  enriched  electrode  materials 
with  a  Mossbauer  isotope  (57Fe,  119Sn)  used  as  a  local  probe  in 
Li4Ti50i2  [7-9],  L^T^Oy  [10,11]  for  instance.  In  other  cases  spec¬ 
tra  are  broadened  due  to  a  mixture  of  phases  or  an  amorphization 
of  the  particles.  The  help  of  electronic  structure  calculation  [12,13] 
allows  non-trivial  experimental  attributions  when  mathematical 
ambiguities  due  to  overlapping  peaks  was  observed.  Interpreta¬ 
tion  of  both  isomer  shift  and  quadrupole  splitting  variations  was 
explained  with  such  theoretical  approaches  [14,15].  In  the  case 
of  glassy  compounds  quadrupole  splitting  distributions  have  been 
used  [16,17]  based  on  previous  works  [18-20]. 

Rechargeable  lithium  batteries  have  proven  to  be  one  of  the 
most  successful  solutions  to  achieve  good  electrochemical  perfor¬ 
mances  in  terms  of  cycleability,  low  cost,  environmental  benefit, 
easy  to  produce  and  safety  in  handling  and  operation.  Phosphate 
materials  (Li3Fe2(P04)3  [21])  have  been  investigated  in  the  past 
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for  their  good  ionics  properties  as  solid  electrolytes.  The  LiFeP04 
compound  has  been  chosen  as  cathode  material  because  of  the 
requirements  [22,23].  Pure  LiFeP04  (triphylite)  and  FeP04  (het¬ 
erosite,  berlinite  [24,25],  high  pressure  [26]  and  amorphous  [27] 
polymorphs)  have  been  earlier  studied  for  their  electrochemi¬ 
cal  performances  [28-30]  and  their  physico-chemical  properties 
like  antiferromagnetism  [31  ]  by  X-ray  [32,33]  and  neutron  [34,35] 
diffraction  techniques,  Mossbauer  spectroscopy  [36-38],  Raman 
[39,40]  and  infra-red  [41-43]  spectroscopies,  X-ray  absorption 
[44]  and  X-ray  photoelectron  spectroscopies  [45,46].  The  origin 
of  capacity  loss  during  the  first  cycle  was  also  discussed  for  this 
good  candidate  as  cathode  material  [47].  Due  to  low  ionic  and  elec¬ 
tronic  conductivity  for  LiFeP04,  and  to  its  ageing  in  presence  of 
moisture  [48,49],  optimized  particle  size  [50],  mineral  coating  or 
substitutions  (LiFei_xMnxP04  [51-54]  or  Lii_xFei+xPi_xSix04  [55]), 
doping  by  transition  metals  [56],  or  carbon  composites  [57]  have 
been  proposed  to  overcome  these  issues. 

In  this  paper  we  illustrate,  taking  advantage  of  2D-correlation 
spectroscopy  analysis,  a  way  to  analyse  Mossbauer  spectra 
recorded  in  operando  conditions.  We  have  used  Python  [58]  pro¬ 
gramming  language  with  Numpy  [59]  and  Pylab  [60]  libraries  to 
get  2D-correlation  spectra.  This  concept  of  2D-correlation  was  ini¬ 
tially  proposed  by  Noda  [61,62].  In  general  2D-correlation  spectra 
are  obtained  using  a  perturbation  (temperature  [63,64],  pH  [65], 
lithium  amount  [66],  kinetic  data  [67])  applied  to  the  pristine  mate¬ 
rial  investigated  by  a  given  or  several  experimental  spectroscopies. 
More  details  on  this  analytical  method  can  be  found  here  [68].  From 
our  knowledge  2D-correlation  spectroscopy  analysis  (2D-COSA) 
used  for  various  experimental  techniques  (infra-red  spectroscopy 
[69,70],  mass  spectral  data  [71],  resonance  light  scattering  [72], 
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X-ray  photoelectron  spectroscopy  [  73  ] )  have  not  yet  used  for  Moss- 
bauer  spectroscopy.  We  present  here  an  application  of  2D-COSA  for 
57Fe  in  situ  Mossbauer  spectroscopy  during  electrochemical  extrac¬ 
tion  of  lithium  from  LiFe0.75Mn0.25PO4  with  olivine  structure. 

2.  Theoretical  basis 

2. 1.  Generalized  2D-corr elation  spectroscopy 

A  set  of  m  experimental  Mossbauer  spectra  are  recorded 
as  a  function  of  a  perturbation  variable  x  (e.g.,  Li  amount 
extracted/inserted  during  electrochemical  tests)  ranging  from 
xmin  to  Xmax-  The  Li  variation  between  two  collected  spectra 
is  Ax  =  (xmax  -  Xmin )/(m  —  1 ).  This  continuous  variation  induces  a 
broadening  of  the  experimental  spectra.  Each  raw  spectrum  is 
folded  and  normalized  giving  a  set  of  m  perturbation-induced 
dynamic  spectra,  y(i/,x)  is  then  transformed  into  a  set  of  20- 
correlation  spectra  by  a  form  of  cross-correlation  analysis.  We  first 
define  the  reference  spectrum  y(v)  of  the  system  to  be  the  averaged 
spectrum  given  by: 

1  m 

y(v)=  ~'Yj(v,x)  (i) 

j= i 

Then  we  need  the  dynamic  spectra  obtained  by: 

y(v,x)=y(v,x)-y(  v)  (2) 

The  intensity  of  2D-correlation  spectrum X(iq  ,v2)  represents  the 
quantitative  measure  of  a  comparative  similarity  or  dissimilarity  of 
spectral  intensity  variations  measured  at  two  different  velocities, 
v\  and  v2,  along  with  a  fixed  interval.  In  a  Mossbauer  spectrum  v\ 
corresponds  to  the  position  of  a  Lorentzian  absorption  line.  In  order 
to  simplify  numerical  manipulation,  X(iq  ,v2 )  is  treated  as  a  complex 
number: 

X(vuv2)  =  S(v-i,V2)  +  iA{vuv2)  (3) 

comprising  two  orthogonal  components,  known  respectively  as  the 
synchronous  S(v  1  ,v2 )  and  the  asynchronous  A(iq  ,v2 )  2D-correlation 
intensities. 

Using  Python  programming  language,  S(v \,v2)  is  calculated  as 
follows: 

S  =  dot(y_dyn.conj().transpose(),  y_dyn)/f loat(m  —  l) 

The  intensity  of  the  synchronous  2D-correlation  spectrum  rep¬ 
resents  the  simultaneous  changes  of  two  separate  spectral  intensity 
variations  measured  at  iq  and  v2  within  the  interval  between  xmin 
and  xmax.  The  autocorrelation  spectrum  a(v j),  defined  by  the  diago¬ 
nal  of  S  matrix,  presents  peaks  called  autopeaks,  whose  magnitude 
is  always  positive.  It  represents  the  overall  extent  of  spectral  inten¬ 
sity  variation  at  the  specific  velocity  vx  within  the  interval  between 
xmin  and  xmax.  Cross-peaks  located  at  the  off-diagonal  positions 
of  a  synchronous  2D-spectrum  represent  simultaneous  changes 
of  spectral  intensities  observed  at  two  different  velocities  v\  and 
v2.  Such  a  synchronized  change  suggests  the  possible  existence  of 
a  coupled  origin  of  the  spectral  intensity  variations.  The  sign  of 
synchronous  cross-peak  becomes  positive  if  the  spectral  intensi¬ 
ties  of  the  two  velocities  corresponding  to  the  coordinates  of  the 
cross-peak  are  either  increasing  or  decreasing  together  within  the 
observation  interval.  However,  the  negative  sign  of  cross-peaks 
indicates  that  one  of  the  spectral  intensities  in  increasing  while 
the  other  is  decreasing. 

The  asynchronous  2D-correlation  spectrum  A(v^,v2)  is  calcu¬ 
lated  in  Python  as  follows: 

A  =  dot(y_dyn.conj().transpose(),  dot(Njk,  spectres_dyn)) 

/f  loat(m  —  l) 


In  this  previous  expression  Njk  is  the  Hilbert-Noda  matrix  cal¬ 
culated  as: 

Njk  =  r-! — -  or  0  if  j  =  k 

The  intensity  of  an  asynchronous  2D-correlation  spectrum 
represents  sequential,  but  not  coincidental,  changes  of  spectral 
intensities  measured  separately  at  V\  and  v2.  Cross-peaks  develop 
only  if  the  intensity  varies  out  of  phase  with  each  other.  The  sign 
of  asynchronous  cross-peak  provide  useful  information  on  sequen¬ 
tial  order  of  events  observed  by  the  Mossbauer  spectroscopy  along 
the  lithium  content.  If  the  signs  of  synchronous  and  asynchronous 
cross-peaks  are  the  same,  the  intensity  change  at  iq  occurs  before 
v2.  If  the  signs  of  synchronous  and  asynchronous  cross-peaks  are 
different,  the  intensity  change  at  iq  occurs  after  v2. 

3.  Experimental 

3.1.  Synthesis  procedure 

LiFe0.75Mn0.25PO4  pristine  compound  was  obtained  by  solid- 
state  reaction.  The  precursors:  Li2C03  (Sigma-Aldrich,  99%), 
FeC204-2H20  (Labosi,  99.99%),  Mnn(CH3COO)2  (Sigma-Aldrich, 
98%)  and  NH4H2P04  (Acros  Organics,  99.99%)  were  first  ball  milled 
during  1  h  30  min  and  then  thermally  treated  in  a  tube  furnace 
under  argon  flow  for  8  h  at  800  °C. 

3.2.  X-ray  diffraction  pattern 

LiFe0.75Mn0.25PO4  pristine  compound  has  been  characterized  by 
X-ray  power  diffraction  with  a  Philips  0-20  diffractometer  using 
Cu  Ka  radiation  (A  =  1.5418  A)  and  a  nickel  filter. 

Crystallographic  reference  data  for  LiFeP04  (triphylite)  has  been 
taken  from  power  diffraction  file  (JCPDS  #81-1173,  a  =  10.332  A, 
b  =  6.010  A,  c  =  4.692  A).  These  values  and  the  orthorhombic  Pnma 
space  group  have  been  used  in  the  initial  stages  of  refinement  of 
the  cell  parameters  for  our  LiFe0.75Mn0.25PC>4  compound. 

3.3.  57 Fe  Mossbauer  measurements 

Mossbauer  spectra  have  been  recorded  in  the  constant  accelera¬ 
tion  mode  and  in  transmission  geometry  on  a  standard  Mossbauer 
spectrometer  composed  of  devices  from  Ortec  and  WissEl.  A 
57Co(Rh)  source  with  a  nominal  activity  of  10  mCi  have  been  used. 
Each  Mossbauer  spectrum  has  been  collected  every  128  min  with 
a  deintercalation  uptake  Ax  =  0.053  Li  corresponding  to  an  electro¬ 
chemical  rate  of  C/40  (1  Li  extracted  in  40  h).  The  source  has  been 
always  kept  at  RT.  Calibration  of  the  velocity  scale  was  done  using 
absorption  lines  of  a  25  p,m  iron  foil.  Isomer  shifts  are  given  relative 
to  a-Fe. 

3.4.  Electrochemical  measurements 

For  Mossbauer  measurements  in  operando  conditions,  we  used 
a  designed  cell  allowing  to  use  it  in  a  transmission  configuration. 
The  cell  has  been  designed  in  such  a  way  that  gamma  rays  can  go 
through  the  Be  windows  which  are  the  current  collector,  the  sample 
under  study  is  in  contact  with  an  aluminum  foil  avoiding  corrosion 
of  the  Be  windows  at  high  potential,  and  the  lithium  foil  used  as 
reference  electrode.  Olivine  has  been  used  as  cathode  material  with 
carbon  black  for  better  electronic  conductivity  (80:20)  vs.  Li  used 
as  anode  with  a  Whatman  paper  as  separator  using  1  M  LiPF6  in 
EC:PC:DMC.  We  used  a  rate  of  C/40,  corresponding  to  1  Li  extracted 
in  40  h  with  a  potential  window  ranging  from  3.2  V  to  4.3  V. 
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Fig.  1.  X-ray  diffraction  pattern  of  pristine  LiFe0.75Mn0.25PO4.  Cell  parameters  deter¬ 
mined  by  LeBail  method  [74]  using  the  Pnma  space  group  are  a=  10.356(6) A, 
b  =  6.025(4)A,  c  =  4.703(3)A. 

4.  Results  and  discussion 

For  LiFe0.75Mn0.25PO4  sample,  XRD  measurements  suggest  pure 
phase  within  the  detection  limit  of  the  apparatus  as  shown  in  Fig.  1. 
Our  pristine  material  present  XRD  lines  indexed  in  the  Pnma  space 
group  with  a  set  of  cell  parameters  a  =  10.356(6)  A,  b  =  6.025(4)  A 
and  c  =  4.703(3)  A.  Refined  cell  parameters  have  been  obtained  by 
pattern  matching  using  FullProf  program  [74].  These  values  are  in 
agreement  with  those  published  earlier  [75]. 

In  Fig.  2  we  have  represented  the  57Fe  Mossbauer  spectrum  of 
the  pristine  material.  The  main  contribution  (91%)  with  an  iso¬ 
mer  shift  (IS)  of  1.20(1)  mm  s-1  and  a  quadrupole  splitting  (QS) 
of  2.96(1)  mm  s_1  correspond  to  Fe2+  in  the  olivine  structure. 
These  values  are  in  agreement  with  those  found  in  the  literature 
[37].  Concerning  the  Fe3+  contribution  detected  in  our  Moss- 
bauer  spectra,  the  hyperfine  parameters  are  IS  =  0.3 1  (4)  mm  s-1  and 
QS  =  0.85(8)  mm  s-1.  These  values  are  not  very  closed  to  those  of 
Li*FeP04(OH)x  (taverite)  [55]  the  reaction  product  obtained  from 
LiFeP04  ageing  under  moisture.  The  presence  of  such  an  impurity 
is  not  excluded  since  about  9%  of  Fe3+  is  present  and  only  90%  of 


Velocity  (mm/s) 

Fig.  2.  57  Fe  Mossbauer  spectrum  recorded  at  room  temperature.  Experimental 
(open  circle)  and  calculated  spectrum  (solid  line)  are  represented.  A  small  contri¬ 
bution  (in  red)  due  to  Fe3+  impurity  (see  text  for  explanation)  is  also  shown  for 
clarity.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  the  article.) 


Fig.  3.  Electrochemical  curve  recorded  during  lithium  extraction  from 
LiFeo.75Mn0.25P04.  The  numbered  labels  correspond  to  the  file  number  of  the 
recorded  57  Fe  Mossbauer  spectra  in  operando  conditions. 

the  Li  can  be  extracted  by  electrochemical  reaction  as  can  be  seen 
in  Fig.  3.  An  other  explanation  of  the  Fe3+  contribution  may  due  to 
Fe203.  From  Mossbauer  data  recorded  at  various  temperatures,  it 
has  been  shown  that  anisotropy  constant  can  be  determined  [76]. 
Hence  for  a-Fe203  57Fe  Mossbauer  spectra  present  rather  a  doublet 
instead  of  a  sextet  at  room  temperature  when  particles  are  nano¬ 
sized.  This  doublet  presents  hyperfine  parameters  ranging  from 
IS  =  0.32  mm  s1  QS  =  0.98  mm  s1  (D  =  4nm)  to  IS  =  0.38  mm  s1 , 
QS  =  0.44  mm  s-1  (D  =  1 8  nm)  depending  on  the  mean  particle  size. 
So,  in  our  case,  both  IS  and  QS  could  be  some  Fe203  formed  at  high 
temperature.  From  linear  correlation  between  QS  and  nanopar¬ 
ticle  size  we  found  a  mean  particle  size  of  D=7±2nm  for  our 
Fe203  impurity.  Recently,  Rho  et  al.  have  mentioned  [77]  the  for¬ 
mation  of  both  Fe2P  (singlet  +  doublet)  and  FeP  (IS  =  0.6(3)  mm  s-1 
and  QS  =  0.8(6)  mm  s-1 )  as  impurity  in  LiFeP04. 

The  electrochemical  curve  shown  in  Fig.  3  presents  two  main 
regions.  Starting  for  [x  =  1  Li,  V=  3.2  V]  in  LixFe0.75Mn0.25PO4  we  can 
extract  about  Ax  =  0.68  Li  (Mossbauer  spectra  #1-13)  on  a  mean 
voltage  of  about  3.48  V.  This  voltage  is  known  to  be  Fe2+/Fe3+  redox 
potential.  Then  from  [x  =  0.32  Li,  V=  3.88  V]  to  [x  =  0.10,  V=4.3V] 
we  observe  a  mean  voltage  at  about  4.02  V,  corresponding  to 
Mn2+/Mn3+  redox  potential  (Mossbauer  spectra  #14-18).  Assum¬ 
ing  that  0.91  Li  have  been  used  (9%  of  non-electrochemically  active 
impurity  observed  in  the  Mossbauer  spectrum,  Fig.  2),  Ax  =  0.68 
Li  and  Ax  =  0.23  Li  are  used  to  oxidize  Fe2+  to  Fe3+  and  Mn2+  to 
Mn3+,  respectively.  These  values  are  rather  close  to  the  expected 
one:  0.75  x  0.91  =0.683  to  oxidize  Fe2+  and  0.25  x  0.91  =0.227  to 
oxidize  Mn2+,  in  a  rather  good  agreement  with  the  stoichiometric 
ratio  Fe:Mn  of  0.75:0.25. 

We  have  represented  in  Fig.  4  the  57Fe  Mossbauer  spec¬ 
tra  obtained  during  electrochemical  Li  extraction  as  mentioned 
above.  We  clearly  see  the  evolution  of  the  absorption  lines 
upon  the  perturbation,  i.e.,  the  Li  content.  In  spectrum  labeled 
1,  we  retrieve  the  doublet  of  the  olivine  (LiFe2+0.75Mn2+0.25PO4) 
constituted  of  the  two  peaks  located  at  v\  =  -0.28  mms-1  and 
y2  =  2.68  mms-1,  corresponding  to  the  hyperfine  parameters 
IS  =  1 .2  mm  s-1  and  QS  =  2.96  mm  s-1 .  After  extraction  of  Ax  =  0.68 
Li  (spectrum  labeled  13),  we  observe  the  doublet  of  the  lithio- 
philite  (Li0.25Fe3+o.75Mn2+0.25PC)4)  constituted  of  the  two  peaks 
located  at  V\  =-0.22  mm  s-1  and  v2  =  0.98  mms-1,  corresponding 
to  the  hyperfine  parameters  IS  =  0.38  mm  s-1  and  QS  =  1 .20  mm  s-1 . 
At  the  end  of  lithium  extraction  (spectrum  labeled  18),  we  observe 
the  doublet  of  the  heterosite  (Fe3+0.75Mn3+0.25PO4)  constituted  of 
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Perturbation-induced  dynamic  spectra 


Velocity  [mm/s] 


Fig.  4.  57Fe  Mossbauer  spectra  recorded  in  operando  conditions  during  Li  extraction 
from  pristine  LiFeojsMnojsPCU.  Labels  1  and  1 8  correspond  to  LiFeo.75Mn0.25P04  and 
Feo.75Mn0.25P04,  respectively.  Colored  lines  [a-g]  correspond  to  cross-correlation 
peaks  positions  observed  in  the  synchronous  spectra.  The  evolution  of  the  intensity 
for  those  specific  positions  is  discussed  in  the  text.  (For  interpretation  of  the  refer¬ 
ences  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the 
article.) 


the  two  peaks  located  at  v\  = -0.36 mm s-1  and  i/2  =  1.18mms_1, 
corresponding  to  the  hyperfine  parameters  IS  =  0.41  mm  s-1  and 
QS  =  1.54mms_1. 

In  Fig.  5  we  have  plotted  the  averaged  spectrum  which  give  an 
indication  of  the  regions  affected  by  Li  extraction.  We  see  again 
three  main  peaks  located  at  about  [-0.2, 1.1  and  2.6  mm  s-1].  This 
averaged  spectrum  is  needed  to  calculate  both  2D-correlation  syn¬ 
chronous  (Fig.  6)  and  asynchronous  (Fig.  7)  spectra. 

In  the  2D-correlation  synchronous  spectrum  (Fig.  6),  we  see 
autopeaks  with  positive  contribution  (diagonal  represented  in 
Fig.  8)  and  cross-peaks  with  positive  and  negative  intensities 
depending  on  the  respective  evolution  of  the  intensities  of  two 
given  velocities.  Positive  and  negative  cross-peaks  have  the  follow¬ 
ing  coordinates:  ([a,  f],  [d,  e]  and  [c,  g])  and  ([a,  g],  [e,  g]  and  [f,  g]), 
respectively.  In  Table  1,  we  have  reported  the  characteristic  values 
of  these  coordinates. 

The  transmission  intensity  is  slightly  correlated  due  to  the  pres¬ 
ence  of  small  positive  cross-peak  with  [c,  g]  coordinates  in  the 
2D-correlation  synchronous  spectrum.  In  Fig.  7,  cross-peak  located 


Averaged  spectrum 


Velocity  lmm/sj 

Fig.  5.  Averaged  spectrum  determined  from  the  mean  of  the  18  Mossbauer  spec¬ 
tra  recorded  during  extraction  of  Li.  The  main  absorption  peak  located  at  negative 
velocity  is  due  to  the  overlap  between  contributions  of  both  Fe2+  and  Fe3+.  The  peak 
located  near  1  mms-1  correspond  to  Fe3+  contribution  and  the  peak  at  2.6mms_1 
correspond  to  Fe2+  contribution. 


at  [c,  g]  presents  a  null  intensity  in  the  2D-correlation  asynchronous 
spectrum.  In  such  analysis,  positive  (negative)  values  mean  that 
variation  at  a  given  velocity  occurs  before  (after)  the  perturbation. 
In  the  context  of  Mossbauer  spectroscopy,  a  close-to-zero  value 
in  the  2D-correlation  asynchronous  spectrum  suggests  that  varia¬ 
tion  observed  at  two  given  velocities  are  coincidental.  Hence,  these 
two  positions  belong  to  the  same  doublet  of  a  specie  defined  by  an 
isomer  shift,  IS  =  (t/i  +  v2  )/2  and  a  quadrupole  splitting,  QS  =  v2  -  V\ . 
Using  this  kind  of  analysis,  we  propose  an  attribution  given  in 
Table  1.  In  Table  2,  we  have  reported  the  hyperfine  parameters 
deduced  from  the  2D-correlation  analysis  and  those  obtained  from 
a  classical  fitting  procedure.  Our  2D-correlation  becomes  a  new  tool 
for  analyses  of  rather  complicated  Mossbauer  spectra  presenting 
some  overlaps  between  absorption  lines  and  small  variation  from 
a  spectrum  to  an  other.  Another  point  to  emphasize  is  that  2D- 


Synchronous  spectrum 


I 


0.00154921 

0.00103280 

0.00051640 

0.00000000 

-0.00051640 

-0.00103280 

-0.00154921 

-0.00206561 


Fig.  6.  2D-correlation  synchronous  spectrum  obtained  from  experimental  57 Fe 
spectra  recorded  during  electrochemical  lithium  extraction  from  LiFeo.75Mno.25P04. 
Positive  (light  green)  and  negative  (dark  green)  cross-peaks  are  observed.  Positive 
cross-peaks  have  the  following  coordinates:  [a,  f],  [d,  e]  and  [c,  g],  while  negative 
cross-peaks  have  mainly  the  following  coordinates:  [a,  g],  [e,  g]  and  [f,  g].  (For  inter¬ 
pretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  the  article.) 
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Fig.  7.  2D-correlation  asynchronous  spectrum  obtained  from  experimental  57Fe 
spectra  recorded  during  electrochemical  lithium  extraction  from  LiFeojsMnojsPCU. 
Positive  cross-peaks  (light  green)  negative  cross-peaks  (dark  green)  suggest  varia¬ 
tion  at  a  given  velocity  occurs  before  (after)  for  another  velocity.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  the  article.) 


Autocorrelation  spectrum 


Fig.  8.  Autocorrelation  spectrum  determined  from  the  2D-correlation  synchronous 
spectrum.  Features  previously  observed  in  both  synchronous  and  asynchronous 
correspond  to  either  close  to  maxima  in  the  autocorrelation  spectrum  [a,  d-f]  or 
minimum  [c].  This  latter  peak  is  more  broadened  than  those  centered  at  about 
2.7  mm  s-1 .  It  contains  some  meaningful  information  as  can  be  read  in  the  text. 


correlation  applied  to  Mossbauer  spectroscopy  allows  to  extract 
both  isomer  shift  and  quadrupole  splitting  of  the  species  using  the 
whole  set  of  spectra  instead  of  using  data  separately. 

The  autocorrelation  spectrum  plotted  in  Fig.  8  shows  the  posi¬ 
tions  of  these  characteristic  values.  These  sets  of  coordinates  are 
used  to  represent  the  relative  transmission  variation  upon  Li  con¬ 
tent  as  shown  in  Fig.  9. 

Relative  transmission  at  velocity  [g]  continuously  increases 
from  file  #1  to  file  #13  or  file  #14  and  then  still  constant  up 


Table  2 

Comparison  of  isomer  shift  and  quadrupole  splitting  (<5caic  and  Z\caic  in  mms"1) 
deduced  from  2D-correlation  analysis  and  those  obtained  by  least-square  fitting 
procedure  (<5fit  and  Afit  in  mm  s_1 ). 


V\ 

V2 

<5Calc 

^  calc 

<5fit 

^fit 

Fe11 

-0.22(c) 

2.68(g) 

1.23 

2.90 

1.20 

2.96 

Fem! 

-0.02(d) 

1.02(e) 

0.52 

1.04 

0.38 

1.20 

Fem2 

-0.36(a) 

1.18(f) 

0.41 

1.54 

0.41 

1.54 

to  the  end  of  the  reaction  (file  #18).  This  unambiguous  posi¬ 
tion  corresponds  to  the  oxidation  process  of  Fe2+  to  Fe3+  in 
LixFeo.75Mn0.25P04.  From  x  =  1  (file  #1 )  to  x  =  0.34  (between  file  #13 
and  file  #14),  we  can  write  the  redox  process  as  the  following: 

LiFe2+0.75Mn2+o  25PO4  -yLi  -*  Li1_yFe2+o.75-yFe3+yMn2+o.25P04 

(0  <y <  0.75)  (i) 

Now,  focussing  our  attention  on  relative  transmission  at  veloc¬ 
ities  [d-f],  we  can  have  an  idea  of  the  mechanism  involved  during 
the  next  steps  of  the  Li  extraction.  In  fact,  it  corresponds  in  some 
way  to  what  happens  in  the  Fe3+  local  environment.  It  is  clear  that 
any  variation  detected  in  the  Mossbauer  spectra  may  be  due  to  oxi¬ 
dation  of  Mn2+  to  Mn3+.  Quadrupole  splitting  is  known  to  be  very 
sensitive  to  second  neighbors  effect,  and  to  the  charge  of  the  ions.  In 
Fig.  9  we  can  see  that  the  relative  transmission  at  [e]  decreases  from 
file  #1  to  file  #1 3  and  still  constant  until  the  end  of  the  Li  extraction. 
While  the  relative  transmission  at  [d]  decreases  from  file  #1  to  file 
#13  and  then  increases  from  file  #1 4  to  file  #1 8.  These  observations 
suggest  the  following  process: 

Li0.25Fe3+o.75 Mn2+0  25 P04  -  (y-0.75)Li 

-»•  Li1_yFe3+o.75Mn2+1_yMn3Vo.75P04(0.75  <y<  1)  (ii) 

Flence,  the  presence  of  Mn3+  in  the  neighborhood  of  Fe3+  induces 
a  higher  quadrupole  splitting  due  to  the  second  neighbor  shell.  This 
can  be  related  to  change  induced  in  the  local  symmetry  of  the  Fe3+ 


Fig.  9.  Evolution  of  the  relative  transmission  for  various  velocities  [a-g]  vs.  file  num¬ 
ber.  A  decrease  in  the  relative  transmission  corresponds  to  an  increase  of  absorption 
at  a  given  velocity. 


Table  1 

Velocities  of  the  main  features  [a-g]  deduced  from  cross-peaks  analysis  of  both  2D-correlation  synchronous  and  asynchronous  spectra  and  autocorrelation  spectrum.  Some 
attributions  are  given  as  help  for  the  discussion. 


Features 

a 

b 

c 

d 

e 

f 

g 

Velocity  (mms-1) 

-0.52 

-0.36 

-0.22 

-0.02 

1.02 

1.18 

2.68 

Attribution 

i7i(Fem2) 

i7i(Fen) 

^(Fe111! ) 

i72(Femi) 

i72(Fem2) 

v2(FeH) 
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site  due  the  Jahn-Teller  effect  on  Mn3+  (3d4).  We  also  observed  this 
by  in  situ  XRD  measurement  (not  the  purpose  of  this  paper)  and 
our  results  are  rather  in  a  good  agreement  with  those  previously 
published  by  Yamada  et  al.  [57]. 

Concerning  the  Mossbauer  investigations,  we  have  shown  by  in 
situ  and  operando  measurements  that  57 Fe  is  directly  sensitive  to 
oxidation  state  of  the  iron  (Fe2+/Fe3+)  during  the  first  stage  of  the 
Li  extraction  as  expected  in  Fig.  3.  We  have  also  shown  that 57 Fe  is 
indirectly  sensitive  the  second  neighbor  shell  through  variations 
of  Mn  oxidation  state.  The  quadrupole  splitting  is  ranging  from 
1.20  mm  s-1  when  Fe3+  are  surrounded  by  Fe3+  and  Mn2+  sharing 
the  M2  site  of  the  olivine  structure,  to  1.54mms_1  when  Fe3+  are 
surrounded  by  Fe3+  and  Mn3+. 

5.  Conclusion 

2D-correlation  spectroscopy  applied  to  Mossbauer  spec¬ 
troscopy  during  electrochemical  studies  of  electrode  materials  in 
real  conditions  (in  situ  and  operando)  can  be  considered  as  a  new 
tool  for  the  Mossbauer  community.  This  analytical  technique  gives 
valuable  information  on  the  electrochemical  mechanisms  involved 
in  extraction  (insertion)  of  lithium  from  (in)  electrode  materials. 
Compared  to  other  spectroscopies,  Mossbauer  spectroscopy  is  a 
complementary  technique  since  this  electrode  material  is  based 
on  Fe2+/Fe3+  redox  system,  57Fe  and  gives  valuable  information  at 
atomic  scale  in  the  bulk  material,  while  XPS  and  Raman  are  more 
dedicated  to  surface  analysis.  With  the  help  of  2D-correlation  anal¬ 
ysis,  Mossbauer  spectroscopy  can  be  used  in  operando  conditions, 
i.e.,  during  Li  insertion  (extraction)  in  (from)  electrode  materials 
tested  in  real  battery  with  the  presence  of  current  collector  (Be 
windows),  separator  and  Li  foil.  Hence  Mossbauer  effect  measure¬ 
ments  yield  the  evidence  for  the  oxidation  processes  taking  place 
in  the  course  of  lithium  extraction. 

2D-correlation  technique  can  be  used  in  other  interesting 
application  like  catalysis,  in  situ  corrosion  measurements,  f- 
Lamb-Mossbauer  factor  determination  for  instance. 
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